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Abstract. We quantify the effect of the spin-orbit interaction on the Rydberg- 
like series of image state electrons at the (111) and (001) surface of Ir, Pt and Au. 
Using relativistic multiple-scattering methods we find Rashba-like dispersions with 
AE so (K) — jK with values of 7 for n — 1 states in the range 38 - 88 meV A. 
Extending the phase-accumulation model to include spin-orbit scattering we find that 
the splittings vary like l/(n+ a) 3 where a is the quantum defect and that they are 
related to the probability of spin-flip scattering at the surface. The splittings should 
be observable experimentally being larger in magnitude than some exchange-splittings 
that have been resolved by inverse photoemission, and are comparable to linewidths 
from inelastic lifetimes. 



Submitted to: J. Phys.: Condens. Matter 

PACS numbers: 73.20.-r, 79.60.Bm, 71.15.Rf, 71.70.Ej 

1. Introduction 

Image states jT] are a special class of weakly bound surface electronic states in which an 
electron outside a dielectric or conductor polarises the surface and is then attracted to 
the resulting "image charge". Asymptotically the potential varies like V(r) ~ —(Az) -1 
so that a band gap preventing penetration of electrons into the crystal leads to a 
Rydberg-like series of states which in the case of a purely Coulombic image force 
at a planar metal surface arise at energies E n (K) = —[0.85 eV]/n 2 + h 2 K 2 /2m, 
n = 1,2... where K is the electron wave vector parallel to the surface and m the 
electron mass. Deviations from this behaviour reflect the influence and response of the 
surface-dependant electronic and atomic structure, which may therefore be investigated 
by studying image states. Examples of theoretical and experimental work include the 
systematics of image states binding and dispersion on clean surfaces [2] , image state on 
overlayers (HHHE], a ^ stepped metal surfaces [Hj, exchange splitting of image states at 
ferromagnets [Zl IH1 IH] , as well as image states at surface nanostructures jTUl HU E21 UH] • 
In recent years there has also been considerable interest in the dynamics of image state 
electron [HI UHl EH Ej as model electronic excitations at surfaces. 

f To whom correspondance should be addressed (s.crampin@bath.ac.uk) 
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One aspect of the physics of image state electrons that has yet to be addressed is 
the influence of the spin orbit interaction H so = [hj '4m 2 c 2 )er ■ (W x p) |ZH] which has 
recently been found to have a significant effect on other surface state electron levels at 
the surfaces of conductors with high atomic number O 120] • At first sight the spin- 
orbit interaction might be expected to be negligible. The mathematical analogy that 
can be drawn between the Schrodinger equation describing the electrons moving in the 
Coulomb-like image potential and that of s-electrons in the hydrogen atom enable the 
image state wavefunctions to be written as 

*Pn,K, s ( r ) = (l/8)zR n0 (z/A) exp(iK ■ r l{ ) Xs (1) 

where R n i{f) is the normalised radial hydrogenic wavefunction and \s a Pauli spinor: 
Xi = (o) > Xi — (?) ■ Using these wavefunctions to diagonalise the spin-orbit perturbation 
H so in the subspace of degenerate image state levels gives a spin-orbit splitting of 

A£f> = ^ 2K = r 012 meV A i * (2) 
64(47re )n 3 L J n 3 K J 

where a = e 2 /47reo^c is the fine structure constant. This is well below the current 
resolution of inverse photoemission, two-photon photoemission or scanning tunnelling 
spectroscopy. However, it has previously been recognised that a more significant 
contribution to the spin-orbit splitting of "crystal-derived" surface states arises from the 
brief time spent by the electron in the vicinity of the nuclei of the surface atoms, where 
the gradient contribution to the spin-orbit interaction is | VV| ~ Z/r 2 . In this paper we 
report on calculations that we have performed to quantify the magnitude of the spin- 
orbit splitting that arises from the penetration of the image state wave function into the 
crystal at surfaces of Ir, Pt and Au. These are described in sectional In section El we 
describe the modification of the phase accumulation model for image state energetics to 
include the effects of the spin orbit interaction. Finally, we summarise and discuss our 
findings. 



2. Relativistic electronic structure calculations 



To calculate the spin-orbit splitting of image states we use a recently developed code 
that implements relativistic multiple-scattering theory. The theory behind this method 
is essentially that described by Halilov et al. [2T] so we do not reproduce it in detail here. 
The basic idea is that the electronic structure is found from the single-particle Green 
function corresponding to the Dirac Hamiltonian H = ca ■ p + (3mc 2 + V ^H|- Thus 
spin-orbit effects are treated non-perturbatively. Using scattering techniques the Green 
function is determined for the special case of semi-infinite crystals with two-dimensional 
in-plane translational periodicity, treating intralayer scattering within an angular 
momentum representation and interlayer scattering in a plane wave representation. Our 
calculations use 25 and 19 two-dimensional reciprocal lattice vectors to describe the 
interlayer scattering for the (001) and (111) surfaces respectively, and partial waves up 
to £ max = 4 22 . The semi-infinite substrate means that continuum and surface- localised 
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Figure 1. Calculated surface state dispersion curves. The shading indicates the 
presence of bulk or vacuum continuum states, (a) Pt(lll) n = 1 and n = 2 image 
states. The inset demonstrates the splitting of the n = 1 state is linear in K. (b) 
Au(001) n = 1 and n — 2 image states, (c) Au(lll) surface state (note the energy 
scale in this case is with reference to the Fermi level). 



states are clearly distinguished in the wave vector resolved local density of states, found 
from the imaginary part of the Green function. As in the non-relativistic version of 
the code |2H] the electronic structure is found self-consistently using the local density 
approximation to density functional theory. We use the atomic sphere approximation 
for the crystal potential (including dipole contributions), with the potential in the three 
outermost atomic layers allowed to vary in response to the presence of the surface. Since 
the local density approximation does not lead to an image-like surface barrier, and hence 
does not support image states, once self-consistency has been achieved we replace the 
self-consistent barrier with a parameterised model barrier, for which we use the "JJJ" 
potential [21] 

V(z) -S (l-e x ^)/4(z-z ), z<z 

Vb{Z) ~ \ -U/ (1 + Ae-K'-*)) , z > zo. [ ' 

The fitting parameters A, U and Zq were fixed by starting with values quoted by Smith et 
al. j2S], who fitted to first-principles slab calculations, and then adjusted slightly to place 
the n = 1 image state at K = 0, Ei, close to values found experimentally. The procedure 
does not uniquely fix the parameters, but we found that different combinations that gave 
the same value for E± resulted in almost identical image state dispersion curves. Note 
that our results are for the lxl unreconstructed surfaces of the materials studied. In 
several cases the surfaces undergo complex surface reconstructions (e.g. Au(001) and 
Pt(001) adopt a 5 x 20 reconstruction but may be prepared in the lxl structure - see 

TO- 

In figure we illustrate the dispersion curves that we find. For both Pt(lll) and 
Au(001) the spin-orbit interaction can clearly be seen to split the n = 1 image state, 
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Table 1. Calculated image state energies including the spin-orbit interaction at (001)— 
1 x 1 and (111) - 1 x 1 surfaces of Ir, Pt and Au: E n {K) = E n + h 2 K 2 /2m ± (j/2)K. 
E n values in brackets are from experiment. 



surface 


n 


E n (eV) 


m(m e ) 


7 (meV A) 


Ir(lll) 


1 


-0.65 


0.95 


56 ± 1 




2 


-0.18 


1.00 


9± 1 


Ir(001) 


1 


-0.61 


0.94 


38 ±1 




2 


-0.17 


0.99 


6± 1 


Pt(lll) 


1 


-0.69 (-0.65°, -0.78 fc ) 


1.05 


50 ±2 




2 


-0.19 (-0.16 a , -0.20 fc ) 


1.03 


8± 1 


Pt(001) 


1 


-0.60 (-0.60 c ) 


0.96 


47 ±2 




2 


-0.17 


0.99 


9± 1 


Au(lll) 


SS 


-0.50 (-0.41 d ,-0.49 e -0.51^) 


0.23 


800 ± 50 


Au(001) 


1 


-0.66 (-0.69V0.63' 1 ) 


1.05 


88 ±4 




2 


-0.18 


1.05 


20 ±2 



a See Ref. [H] 
b See Ref. 30, 
c See Ref. [2HI 
d See Ref. [E] 



e See Ref. [3T] 
f See Ref. [33] 
9 See Ref. [3J 
h See Ref. [35] 



and whilst a splitting exists for the n = 2 and higher states it is much smaller. The 
inset in the figure ^ illustrates the variation of the splitting with wave vector K , the 
near-linear variation corresponding to a Rashba-like dispersion [271 [2E] 

E n (K) ~E n + ^ ± ( 7 /2)tf. (4) 

From curves such as these we extract Rashba-parameters 7 by a least-square fit using 
wave vectors K < 0.2 A . These values are tabulated in table Qfor the image states at 
the (111) and (001) surfaces of Ir, Pt and Au. At Au(lll) the vacuum level lies outside 
of the projected band gap so that the image states in this case exist as resonances. We 
have not therefore included results for this case, but instead give the results that we find 
for the spin orbit splitting of the occupied surface state that occurs within the band gap 
at this surface. The dispersion of this state is shown also in figure and agrees well 
with previous work, validating our procedure. We find the wave vector splitting at the 
Fermi energy is Ak = 0.023 A , compared with experimental values of 0.023 A [19J, 
0.025 A" 1 and 0.027 A~ X HE], and a theoretical value of 0.023 A" 1 [37| that 

have been reported previously. 

The results in table [T] indicate that the spin-orbit splitting of n = 1 image states at 
Ir, Pt and Au surfaces is an order of magnitude smaller than that of the Au(lll) 
Shockley surface state, and that of the n = 2 states smaller still. These trends 
reflect the differing extents to which the corresponding wave functions penetrate the 
crystal and experience the spin-orbit interaction at the ion cores. At the (001) surfaces 
AE^l < AE^ 1 < AE^ U which might be expected given the increasing atomic number 
(Z Ir = 77,Z Pt = 78,Z Au = 79), but at the (111) surface AE l * > AE%\ pointing to 
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a more complicated "band-structure" effect. In figure Q the splitting of the Au(OOl) 
n — 1 image state is also seen to be affected as it disperses towards the band edge of 
continuum levels. We return to this later. 



3. Phase accumulation model 



The "standard model" for understanding image state energies is the phase accumulation 
model PP in which surface states are envisaged as one-dimensional waves trapped by 
multiple reflection from the surface barrier and the crystal. In a region of constant 
potential (see figure |2J between barrier and crystal (which may be infinitesimal in width) 
the electron wave function can be expressed in terms of forward and backward travelling 
plane waves 

if)(z) = Aexp(ikz) + B exp(— ikz). (5) 

The two components are related at the barrier reference plane z = zb by the barrier 
reflection coefficient r B , 4>(z) oc exp(— ik(z — z B )) +r B exp(iA;(,2 — z B )) and at the crystal 
reference plane z = Zq by the crystal reflection coefficient rc, ip{z) oc exp(ik(z — zq)) + 
r c exp(— ik(z — Zq)) which together with © give rise to the condition for a surface state 
to exist: 

t-qTc exp(2i/cd) — 1 = 0, d = zb — zq. (6) 

For energies below the vacuum level and coincident with the crystal band gap the 
reflection probability at both crystal and barrier are unity and the reflection coefficients 
may be written in terms of phases: r B = exp(i0B), = exp(i0c)- The surface state 
condition then becomes 

0b + 0c + 2kd = 2vrn n = 0, 1, . . . (7) 

which is a condition on the round-trip phase accumulated by the electron wave. The 
phases in flj|) increase with energy. The crystal phase increases from to 7r as the energy 
sweeps across the band gap, changing most rapidly near the band edges. Towards the 
bottom of the gap it is this variation in <ftc which will determine whether or not (JJJ) is 
satisfied, so that any surface state that does arise is usually referred to as "crystal- 
derived". On the other hand, 0b increases more and more rapidly as the energy 
approaches the vacuum energy, varying to a good approximation as 



ME)=^[\r-^-A- (8) 

In combination with (JJJ) this yields a Rydberg-like series of image states 

B „ = Z^X „ =1 , 2 ,... (9) 
(n + a) z 

where the quantum defect a = (1 — [0c + 2kd]/Tc) /2 may usually be considered constant 
over the range of energies at which the image states are found. 
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Figure 2. Schematic illustration of quantities entering the phase accumulation model. 



We now consider the extension of this model to include the spin orbit interaction. 
Introducing the electron spin in to the wavefunction in (JHJ 

^K^D-GO^KS"^- < io > 

and reflection from the crystal is now described by a matrix 

«-(:s:s) (") 

allowing for the possibility of spin-flip upon reflection. With a similar matrix used to 
describe scattering from the barrier, the condition for a surface state becomes 

det [R B R C exp(2ifcd) - 1] = 0. (12) 

The four reflection coefficients in (jllj) are not independent - for example flux 
conservation requires that \rj]\ 2 + \tq\ 2 = 1 within a gap, and for a non-magnetic 
crystal rjl = Tq. In the non- magnetic case and for a planar potential V = V(z) the 
electron wave functions ^K{ r ) — iPk{ z ) exp(iif • rn) are found from the Hamiltonian 

H = -f-V 2 + V+ -^<r ■ ( W x K) . (13) 

which may be diagonalised by a rotation in spin space 

-, h 2 , h 2 K alV , s 

H' = U t HU t - = r + V + —- ns (14) 



with 




_ ■ , 1 -iexp(-M) . 
U »- — \ 1 +iexp(-h?) 1 (15) 
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where $ is the angle of the electron wave vector. The Hamiltonian H' does not mix spin- 
up and spin-down channels and so in this representation the reflection matrix describing 
scattering from the crystal is diagonal: 

Since the spin-orbit interaction is negligible in the barrier the barrier reflection matrix 
is also diagonal, and for a non-magnetic surface = = = = r B so that 
the surface state condition (fT2~j) becomes r^r-Q exp(2ifcd) = 1 leading to the round-trip 
phase condition 

0c ± V + 0b + 2kd = 2nn n = l,2,... (17) 

where we have introduced = exp(i(0c±?7)) appropriate to energies within a gap. The 
surface states now come in spin-split pairs, with spin orientations that may be deduced 
from the spinors that are obtained by rotating back in to the original reference frame 
the spin-up and -down eigenspinors x' + = ( ) > x'~ = (i) °f the primed frame: 

^^^UxW))' <18) 

Thus S = ±(— sin?9, cos$, 0) = ±z x .K" and we find that the spins lie in the surface 
plane and perpendicular to K, to the left (right) for + (— ). 

Rotating back to the original spin frame gives the reflection matrix (jllj) as 

d tt-i&tt ca\( cosr l exp(-h?)sin?7 \ 

R c = U# R C U# = exp i0c . 19 

\ — exp{iv) smr] cosr] I 

In figure El we illustrate the variation in tan?? at the Au(001) surface calculated from 
the reflection matrix found using the relativistic multiple-scattering method of section 
121 For K along FX (see inset in figure 0J) equation f)19|) tan?? = r^^/r*'. The 
relativistic multiple-scattering calculations include the atomic structure of the surface 
and the crystal potential is not one-dimensional, but the reflection coefficient for 
specular reflection behaves in a very similar manner to that of a one-dimensional crystal, 
especially for small K where the wave function varies only slowly across the surface. In 
particular we find that the angular variation predicted by equation (|19J) is satisfied to 
within a percent or so. It is evident from figure El that the magnitude of the spin- 
orbit induced phase change is small, and away from the band edges rj is approximately 
independent of energy at Au(001) and linear in K: rj^-^JEjK) ~ [—0.25 A]K. 

To first order the round-trip phase condition (|17|) is satisfied at energies E n ± 
AE^°/2 where 

) , (20) 

neglecting the energy dependence of r\ which is small compared to the other phases. 
For all the surfaces that we have studies we have found that r] does not change sign 
across the band gap, and the denominator in (|2T)|) is positive. Hence the model predicts 



AAt () - ^(0B + 0c + 2foZ) 
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Figure 3. Variation of the spin-flip to non-spin-flip scattering ratio r"^/r'' for 
Au(OOl) in the major band gap of the projected band structure, (a) as a function 
of energy at K — (0.2,0.0) A . (b) as a function of wave vector along n at the 
mid-gap energy. The inset shows the projected gaps across the Brillouin zone at the 
same energy. 
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Figure 4. Spin orientation of spin-orbit split surface states for i] < 0. 



a series of spin-orbit split states with identical spin orderings, which we have confirmed 
is also the case in the relativistic multiple scattering calculations described in section 
121 In particular with r\ < the surface states that exist are split with the lower of each 
pair of levels having the spin pointing to the right of K, as shown in figure 0] This is in 
agreement with the spin assignments shown in Henk et al. jSZ] for the Au(lll) surface 
state, but disagrees with those given in Ref. pH] . 

When the gap contains image states, over the narrow range of energies within which 
the image states are found the energy- dependence of the round-trip phase is dominated 
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by the variation in the barrier phase (JSJ) and then 




r] x 1.7 eV 



7r(n + a) 3 



(21) 



Thus the spin-splittings exhibit the same scaling as the lifetime broadening [T, in each 
case the behaviour ultimately originating in the variation with n of the wave function 
overlap with the crystal. We also see from (|2"TJ) that the linear-in-.fr behaviour of AE^° 
arises from similar behaviour in rj. Since r\ is small tan 77 ~ sin 77 ~ 77 to a good 
approximation, and hence 77 ~ — |r^|. Thus the spin-orbit splitting is directly related 
to the spin-flip scattering rate, which could therefore be determined from experimental 
values of image state splittings. At Au(OOl), E\ ~ —0.66 meV (tabled)), so combining 
equations (JSJ) and (JTTJl gives 



and a splitting of 22 meV (figure^) at K = 0.25 A -1 yields rj ~ 0.06, which agrees with 
the value found from the multiple-scattering calculations shown in figure EJ Finally, we 
note that equation (}2*T]) will not hold near band edges where the energy-dependence of 
the crystal phase 0c cannot be neglected. This is the origin of the anomalous dispersion 
shown for the n = 1 state at Au(001) in figure ^ 

4. Discussion 

To summarise, we have investigated the effect of the spin-orbit interaction on image state 
electrons at the (111) and (001) surfaces of Ir, Pt, and Au. Non-perturbative calculations 
that use relativistic multiple-scattering theory with self-consistent potentials and a 
parameterised surface barrier predict Rashba-like dispersion of the image state bands 
with splittings for n = 1 that are a factor 10-20 times smaller than that of the Au(lll) 
Shockley state, for which our results are in good agreement with experiment and previous 
theory. Extending the phase accumulation model to include spin-orbit scattering, we 
find that the splittings scale as 1/(77, + a) 3 , where a is the quantum defect, and are 
directly related to the spin-flip scattering rate at the surface. The largest image state 
splittings we find are at Au(001) for which AEf° = 22 meV at K = 0.25 A -1 . This is 
larger than some exchange splittings of image states that have previously been resolved 
(e.g. 18 ± 13 meV at Ni(lll) and 13 ± 13 meV at Ni(001) 38 ) by exploiting 
the spin resolution available in spin-resolved inverse photoemission, suggesting that the 
spin-orbit splitting may also be observable with an appropriate experimental set-up. 
The splittings are comparable to lifetime broadenings of late-transition and noble metal 
image states [THEE], indicating that account of spin-orbit effects may be necessary when 
determining image state lifetimes from lineshape analysis for 5d metals. 
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